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Abstract 

We propose a modification of tlie standard van der Pauw metliod 
for determining the resistivity and Hall coefficient of flat thin samples 
of arbitrary shape. Considering a different choice of resistance mea- 
surements we derive a new formula which can be numerically solved 
(with respect to sheet resistance) by the Banach fixed point method for 
any values of experimental data. The convergence is especially fast in 
the case of almost symmetric van der Pauw configurations (e.g., clover 
shaped samples). 
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1 Introduction 



The van der Pauw four probe method is a standard technique for measuring 
the resistivity of flat thin samples of arbitrary shape [H |2] . The sample have 
to be homogeneous, isotropic, of uniform thickness and simply connected (i.e., 
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without isolated holes). Four contacts placed on the sample are required. 
They have to be geometric points located on the boundary of the sample (or, 
in practice, errors caused by their finite size should be sufficiently small). 

The van der Pauw geometry is very popular in electric measurements 
and found a lot of applications in physics, compare, e.g., [SI IH |5l El [7] . The 
method consists in performing direct measurement of resistances -Ri2,34 and 
-^23,14 (for more details see the next section), and then using the formula 

exp^ ^J+^"H 

for computing the resistivity p and sheet resistance Rs = p/d of the sample 
of thickness d. Then, the Hall coefficient is computed as 

where A_R24,i3 is the change of -R24,i3 due to the magnetic field B. Equation 
([T]) is believed to be unsolvable by the fixed point method. Usually, instead 
of numerical procedures, a graph of the so called geometric factor is used 
to determine a solution of ([1]). Some authors recommend to use tables of 
numerical values of this function [8] . An inherent inaccuracy of these methods 
seems to be commonly recognized. 

Many attempts have been made to develop and improve the van der Pauw 
approach, see P[Tni[IIl[12l[13l[Tl[l5l[l6]. However, the formula has 
always been treated as a starting point. In this Letter we will show that 
another formula, namely: 

irdRmax f T'd\R2i,13\\ 1 

exp exp = 1 



P \ P / (3) 

(where Rmax = max{i?i2,34, i?23,4i}) , 

can be used instead of ([T]). We will show that preconditions for the Banach 
fixed point theorem are rigorously satisfied for any set of experimental results, 
usually with an excellent rate of convergence. 

Our approach is especially convenient in Hall effect measurements with 
symmetric (or almost symmetric) van der Pauw configuration (e.g., in the 
shape of a clover leaf). Then -R24,i3 is much smaller than Rmax and we need 
few iterations to get very accurate numerical results. 
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2 A brief review of the van der Pauw method 

The main idea of the van der Pauw approach is simple and beautiful. First, 
one considers a sample in the form of the complex upper half plane (with con- 
tacts placed on the real axis). All computations can be explicitly done in this 
case. Then, one applies a deep mathematical theory (the Riemann mapping 
theorem) showing that any other (simply connected) sample is conformally 
equivalent to the upper half plane pTj. What is more, this conformal trans- 
formation preserves all equipotential lines, current lines and boundary condi- 
tions [1]. Therefore any formula which does not contain explicit information 
about positions of contacts is invariant with respect to such transformations, 
and results obtained in the case of the half plane are exactly valid for samples 
of arbitrary shape (provided that they have no isolated holes). 

Therefore, we consider the upper half plane, parameterized by complex 
coordinate z (Imz ^ 0). Four contacts are represented by xi,X2,X3,X4 lying 
on the real axis. In order to perform a measurement we inject electric current 
Jj-fc at contact Xj, take it out at {k ^ j), and measure the voltage between 
remaining two points. Elementary considerations (based on the superposition 
principle) show that electric potential at z is given by 



(note that \zi — Z2\ is a distance between complex numbers zi and Z2). There 
are 4!=24 different ways to perform measurements described above. In any 
case we compute a resistance 



where j, k,m,n are pairwise different (a permutation of 1, 2, 3, 4) and it is 
convenient to denote Rg = p/d {sheet resistance). Thus we have 24 relations 
between xi,X2,X3,Xi and Rg, treated as unknowns. Rjk,mn are calculated 
directly from experimental data. Eliminating xi,X2,X3,X4 van der Pauw 
obtained equation ([1]) valid for samples of arbitrary shape, compare [U [2] . 
We stress that the exact placement of contacts on the circumference of the 
sample is not important with exception of their ordering. 

In the next section we study consequences of equations in more detail. 
In particular, we derive new equation 




(4) 





(5) 
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3 Modification of the van der Pauw method 

In the formula ([5]) one can recognize the cross ratio, a well known and very 
important notion in projective geometry. The cross ratio of four (ordered) 
points Xj,Xk,Xm,Xn is defined as 

The same formula applied for a 4-tuple of complex numbers is used in con- 
formal (Mobius) geometry [TTlllH]. There exists a natural generalization of 
the cross ratio on points in Euclidean spaces of any dimension [T9] . 



Taking into account ([6]) we rewrite equation ([5]) as 

'^Rjk,mn Rs 111 I (-^j; -^fc) •^mi ■^n) \ • (7) 

Cross ratios corresponding to various permutations of four points Xi,X2,X3, 0:4 
are related by a set of identities which can be shortly written as: 

, Xfc, Xni) 2^n) (-^m) '^Jii -^j ; '^/c) (-^j ; -^fci -^n; -^m) ; (8) 
(Xj , Xfc, ~l~ {^■^ j 1 ■^m} -^ki ■^n) 1 ; (9) 

(they can be verified by straightforward elementary calculation). In particu- 
lar, on use of ([8]) and ([9]) we easily derive the following equations 

(Xi, X2] X3, X4)~^ + (X2, X3; X4, Xi)~^ = 1 , (10) 

(Xi, X2; X3, X4) + (X2, X4; Xi, X3) = 1 , (11) 
(X2, X3; X4, Xi) + (X2, X4; Xi, Xs)"^ = 1 , (12) 

Taking into account (I7l), and assuming (without loss of the generality) 

xi < X2 < X3 < X4 , (13) 

we obtain corresponding identities for resistances Rjk,mn- Equations ([8]) yield 
the so called reciprocal and reversed polarity identities, for instance: 

-Rl2,34 = -^34,12 = -^21,43 = -^43,21 • (14) 
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They are useful for eliminating some side effects (one takes an average of the 
above four measurements insted of -Ri2,34, etc.). In our approach improve- 
ments of this kind can be done in exactly the same way as in the standard van 
der Pauw method. Note that ordering f lT3|) means that contacts xi,X2,X3, 
are placed in exactly this order (counterclockwise) on the circumference of 
the sample. 

Cross ratios are not necessarily positive. Using ([H]) and f ll3p we can de- 
termine signs of cross ratios. Moreover, equation ( ITO|) implies upper bounds 
on both (positive) components. Thus: 

(xi,X2;X3,X4) > 1 , 

(x2,X3;x4,a;i) > 1 , (15) 
(a;2,X4;xi,a;3) < . 

Equation (fTOj) yields van der Pauw's formula ([1]). Surprisingly enough, equa- 
tions (ITT]) . ( IT2|) lead to new, physically meaningful, formulas: 

exp(7ri?i2,34/i?s) -exp(7ri?24,i3/^s) = 1- (16) 



exp(7ri?23,4i/-Rs) - exp(-7ri?24,i3/-Rs) = 1- (17) 

For further analysis we choose the first equation if -R24,i3 > or the second 
equation if -R24,i3 < 0. In the first case we have -Ri2,34 > -R23,4i > 0, while in 
the second case -R23,4i > -Ri2,34 > 0. Both cases can be shortly represented 
as equation (E]) where Rmax denotes greater of two values: i?i2,34 or -R23,4i- 



4 Fast converging numerical iterations 

Equation ^ can be rewritten as: 

x = hi{l + e^^) , fc = i|M^, (18) 

where x = nRmax/Rs- The discussion at the end of the previous section 
shows that ^ k < 1. 

Equation f|T8|) has a form x = F{x), characteristic for the Banach fixed 
point method. In order to obtain a solution (the fixed point of the map 
F) one has to iterate: Xn+i = F{xn)- We are going to show that function 
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F{x) = In (l + e^^) satisfies preconditions for tfie Banacli fixed point theorem 
(for any k). Indeed, F maps segment = [in 2, ^^1] into itself because: 

x^ln2 F(x) ^ ln(l + 2'=) ^ ln2 , 

ln2 , , A;ln2 , ln2 (^9) 

wliere we took into account F{x) = kx + ln(l + e~^^'). Tfien, 
k 

l + ~e 

for any x G L^. Therefore, by virtue of Lagrange's mean value theorem 

|F(xi)-F(x2)| 



F'i^)\ = V-ZZj^^k (20) 



Fl — X2\ 



\F'{c)\^k<l (21) 



(for any Xi,X2 G L^) which means that F is a contraction of the segment L^. 

In order to estimate the number of iterations N needed to obtain a pre- 
scribed accuracy 6 we require that the length of the segment after applying 

contractions is smaller than 6: 



The actual number of iterations is, of course, much smaller. Table [1] shows 
the number of iterations needed to obtain the accuracy 6 = 10~^. For k 
approaching 1 the number of iterations increases (tending to infinity). In 
this region {k ^ 1) it is better to use another iterating scheme, see below. 
Note that as an initial point we took xq = In 2 (this is almost obligatory 
for small k, when the length of segment is very small and only a;o = In 2 
belongs to any Lk). Tabled] contains also corresponding values of the relative 
sheet resistance Rg defined by 

^'=R^=x^ (23) 

^ ^m.n.T. ^ 



where x is the solution of ( ITSl) . 

Multiplying equation = 1 + e^^ (equivalent to ( |T8l) ) by we get: 
= 1 — e^^~^. Hence we have another form of equation (1181) : 

x = - ln(l - e-'^'^) , k' = l-k . (24) 
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Table 1: Number of iterations necessary to obtain solution x of Eq. f ITS]) 
(xo = ln2, 6 = 10^5) and Rg = Rg/Rmax as a function of = | -R24, 13 1 /-Rmaa:- 



k 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


N 


1 


4 


5 


6 


7 


9 


12 


15 


24 


42 


Rs 


4.532 


4.302 


4.062 


3.811 


3.546 


3.264 


2.960 


2.623 


2.234 


1.743 


Table 2: Number of iterations necessary to obtain solution x of Eq. (124|) 
(xq = — In /c', 5 = 10^5) and Rg = Rg/Rmax as a function oi k' = 1 — k. 


k' 


0.2 


0.1 


0.01 


10-3 


10-^ 


10-6 


10-8 


10-10 


10-12 


10-15 


N 


22 


18 


10 


8 


7 


5 


4 


4 


4 


3 


Rg 


2.234 


1.743 


0.924 


0.598 


0.434 


0.276 


0.200 


0.157 


0.129 


0.101 



One can rigorously show that preconditions for the Banach fixed point method 
are satisfied (at least for sufficiently small k', namely k' < 0.125) provided 
that as a starting point we take Xq = — In k' (in practice, the Banach method 
seems to work very well for larger range of k', at least up to k' ^ 0.25). We 
omit technical details. Instead, we present Table |2] showing that for small k' 
(i.e., k ~ 1) equation ( 12^ is excellently solvable by the fixed point method. 



5 Summary 

In this Letter we proposed an alternative approach to the standard van der Pauw 
method. Measurements are essentially the same as in the standard method and 
produce three resistances: i?i2,34, -^23,41, -^24,13 (reciprocal and reversed resis- 
tances can be used for improving the accuracy, compare ([HI)). We take i?24,i3 and 
greater of remaining two resistances, denoting it by Rmax- Then we compute two 
coefficients: k = I-R24, 13 1 /-Rmax and k' = 1 — k. In order to find the sheet resistance 
we solve either ([18]) (for ^ /c < 0.9) or ([2lD (for 0.8 < A; < 1) and calculate 
Rg = irRmax/x. In the indicated ranges of k both equations are solvable by the 
Banach fixed point method with excellent rates of convergence. 

Acknowledgments. I am grateful to Kamil Lapihski for turning my attention 
on the van der Pauw method and to Krzysztof Szymahski for discussions. 



7 



References 

[1] L. J. van dcr Pauw: A method of measuring specific resistivity and Hall effect of 

discs of arbitrary sliape, Philips Research Reports 13 (1958) 1-9. 
[2] L. J. van der Pauw: A method of measuring the resistivity and Hall coefficient on 

lamellae of arbitrary shape, Philips Technical Review 20 (1958) 220-224. 
[3] P. Blood, J. W. Orton: The electrical characterisation of semiconductors. Reports 

on Progress in Physics 41 (1978) 157-258. 
[4] Z. Moron: Investigations of van dor Pauw method applied for measuring electri- 
cal conductivity of electrolyte solutions: Measurement of electrolytic conductivity. 
Measurement 33 (3) (2003) 281-290. 
[5] O. Bierwagen, R. Pomracnkc, S. Eilers, W. T. Masselink: Mobility and carrier density 
in materials with anisotropic conductivity revealed by van der Pauw measurements, 
Physical Review B 70 (2004) 165307. 
[6] Y.-W. Tan, Y. Zhang, K. Bolotin, Y. Zhao, S. Adam, E. H. Hwang, S. Das Sarma, 
H. L. Stormcr, P. Kim: Measurement of scattering rate and minimum conductivity 
in graphene, Physical Review Letters 99 (2007) 246803 (4pp). 
[7] J. R. Williams, C. M. Marcus: Snake states along graphene p-n junctions. Physical 

Review Letters 107 (2011) 046602 (4pp). 
[8] A. A. Ramadan, R.D. Gould, A. Ashour: On the Van der Pauw method of resistivity 

measurements, Thin Solid Films 239 (2) (1994) 272-275. 
[9] W. Versnel: Analysis of symmetrical Van dcr Pauw structures with finite contacts, 
Solid-State Electronics, 21 (10) (1978) 1261-1268. 
[10] S. Hyun, M. F. Thorpe, M. D. Jaeger, B. Golding: Resistivity determination from 

small crystallites. Physical ReviewB 57 (11) (1998) 6697-6705. 
[11] W. K. Chan: On the calculation of the geometric factor in a van der Pauw sheet 
resistance measurement. Review of Scientific Instruments 71 (10) (2000) 3964-3965. 
[12] J. D. Weiss, R. ,J. Kaplar, K. E. Kambour: A derivation of the van der Pauw formula 

from electrostatics, Solid-State Electronics 52 (1) (2008) 91-98. 
[13] O. Kouba, A. Nader: Explicit form of the geometric factor for a thin disks resistivity 

measurement. Measurement 43 (9) (2010) 1192-1196. 
[14] J. D. Weiss: Generalization of the van der Pauw relationship derived from electro- 
statics, SoUd-State Electronics 62 (1) (2011) 123-127. 
[15] J. Nahlik, I. Kasparkova, P. Fitl: Study of quantitative infiuence of sample defects on 
measurements of resistivity of thin films using van der Pauw method. Measurement 
44 (10) (2011) 1968-1979. 
[16] H. Li, Y. Sun, W. Wang, H. Hutchinson: Neurocomputing van der Pauw function 
for the measurement of a semiconductor's resistivity without use of the learning rate 
of the weight vector regulation. Journal of Semiconductors 32 (2011) 122002 (8 pp). 
[17] L. V. Ahlfors: Complex Analysis, McGraw-Hill, New York 1953. 
[18] U. Hertrich-.Jeromin: Introduction to Mobius differential geometry, Cambrigde Uni- 
veristy Press 2003. 

[19] J.Cieslinski: "The cross ratio and Cliff'ord algebras", Advances in Applied Clifford 
Algebras 7 (1997) 133-139. 



8 



